begin
  fi = addfile("../power_uv_mcv.nc", "r")
  div = fi->div
  vor = fi->vor

  wks = gsn_open_wks("ps", "../div_vor_global")

  res = True
  res@cnFillOn = True
  res@cnLinesOn = False
  res@mpCenterLonF  = 180

  plot = gsn_csm_contour_map(wks, div, res)
  plot = gsn_csm_contour_map(wks, vor, res)

end
